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Anisotropic eddy viscosity models 

By D. Carati 1 AND W. Cabot 2 


A general discussion on the structure of the eddy viscosity tensor in anisotropic flows 
is presented. The systematic use of tensor symmetries and flow symmetries is shown 
to reduce drastically the number of independent parameters needed to describe 
the rank 4 eddy viscosity tensor. The possibility of using Onsager symmetries for 
simplifying further the eddy viscosity is discussed explicitly for the axisymmetric 
geometry. 


1. Introduction 

Contrary to most of the works presented in this volume, this note does not re- 
sult from a planned project for the summer program. It developed instead from 
discussions during the course of the workshop by many participants concerning the 
representation of anisotropy in the modeling of the subgrid-scale stress in Large 
Eddy Simulation (LES). This study is thus an attempt to present a systematic dis- 
cussion of the influence of anisotropy on the structure of the eddy viscosity tensor. 
Some of the results presented here are not really original since they have been de- 
rived in other contexts (viscoelastic media or magnetized plasmas). However, we 
found several motivations for reproducing the general study of tensor symmetries 
in the special case of the eddy viscosity tensor. 

First, we remark that there is often evidence of anisotropy at the subgrid level. 
The most obvious case arises when the grid itself is anisotropic. In that case, even 
if the flow does satisfy the classical local isotropy assumption, the subgrid velocity 
would be anisotropic by construction. Since most LESs use a non-uniform grid 
with anisotropic stretching, the effects of anisotropy should be taken into account 
in a very wide class of problems. 

Second, the discussions we had during the workshop showed that few attempts 
have been made to introduce the anisotropy at the tensor level in the relation be- 
tween the subgrid scale stress and the resolved strain tensor. On the contrary, most 
of the studies on the influence of anisotropy have focused on possible modifications 
to the isotropic eddy viscosity amplitude (Deardorff, 1970, 1971; Scotti et al . , 1993). 

Finally, the development of the dynamic procedure (Germano, 1992; Ghosal et 
a/., 1995; Lilly, 1992) allows the introduction of multi-parameter models for the 
subgrid scale stress. Therefore, there is no practical reason for practitioners to limit 
their models to an isotropic eddy viscosity. 
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2. Anisotropic eddy viscosity 

In this work we only consider the subgrid scale modeling of an incompressible 
fluid. If the exact description of the large scale pressure is not required, the trace 
of the subgrid scale tensor may be added to the pressure, which is then calculated 
in order to ensure the incompressibility. The only tensor that needs to be modeled 
is 


7 ~*j = u i u j ~ u,u } - - (u k u k - u k u k )6,j . (2.1) 

The usual modeling procedure consists in giving an expression for r,’ in terms of the 
spatial derivatives of the resolved velocity field d t u } . These quantities are usually 
decomposed into a symmetric resolved strain tensor, 

Sij = \ (diUj + d } %ii) , (2.2) 

and an antisymmetric resolved rotation tensor, 

R tJ = ^ ( d,Uj - djtii) = ^€ ijk w k , (2.3) 

where w k is the vorticity and €,j k is the Levi-Civita fully antisymmetric tensor with 
fi23 = 1. The most general tensorial relation in an anisotropic system thus reads: 

T ij — VijkiSki + VijittRki. (2.4) 

For three dimensional turbulence, a naive analysis of this relation would lead to the 
conclusion that both v and (i are described by 81 independent parameters. However, 
very_strong simplifications can be derived by using the tensor symmetry properties of 
T iji Sij an d Rij, as well as the symmetries of the flow. These simplifications do not 
require any assumption (as far as the model (2.4) is accepted). A more debatable 
simplification might apply if the Onsager reciprocal symmetries (Onsager, 1931) are 
assumed to hold for the eddy viscosity tensors. This will be discussed at the end of 
this section. 

2.1 Tensor symmetries 

The tensors r,* and S t] are symmetric and traceless while the tensor R tJ is an- 
tisymmetric. This implies that the eddy viscosity tensor v l]k i has the following 
properties: 


= t'jikl , 
Vijkl = Vijlk , 
Viikl - 0 , 
Vijkk = 0 . 


(2.5) 


Thus, for a given value of (k,l) = ( k*,l * ), the matrix = v ijk .,. is traceless and 
symmetric. Consequently, it has 5 independent components. Similarly, for a given 
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value of(i,j) = the matrix 6*/ = i/^j* is also traceless and symmetric. The 

full tensor u tJ ki is thus described by 5 x 5 = 25 independent parameters. The same 
analysis can be performed for the tensor fiijkh which has the following symmetries: 

f^ijkl “ j^jikl i 

f^ijkl ~ f^ijlk > (2-6) 

j^iikl = 0 . 

Now, the tensor piijkl is symmetric and traceless for its first two indices, while 
it is antisymmetric for its last two indices. Consequently, the full tensor fi t jki is 
described by 5 x 3 = 15 independent parameters. 

2.2 Flow symmetries 

This 25+15 parameter eddy viscosity tensor may be strongly simplified by using 
the symmetries of the flow. Let us consider some simple cases. 

2.2.1 Isotropic turbulence 

Any isotropic tensor can only be constructed with the unit tensor 6ij. Thus, the 
most general isotropic tensor of rank 4 can be written as follows: 


Tijki = ai6tj$ki + 0 , 2 $ik8ji + ozSuSjk . (2.7) 

If we impose the symmetry relations (2.5), it turns out that the eddy viscosity 
tensor v reduces to 


^ijkl — O 


(^>ik$jl + fiiltijk 



( 2 . 8 ) 


while the symmetry relations (2.6) imply that the tensor /j vanishes. Consequently, 
the subgrid scale stress reads: 


r* = -2aSy, (2.9) 

where a is the usual isotropic eddy viscosity (Smagorinsky, 1963). 

The simplest anisotropic situation arises when only one direction can be distin- 
guished from the other. This axisymmetric geometry is thus characterized by a 
vector pointing to the anisotropy direction. We will show that the nature of this 
vector will strongly affect the structure of the eddy viscosity tensor. In particular, 
anisotropy induced by a pseudovector (like a magnetic field or a rotation) must 
be treated differently from the anisotropy induced by an axial vector (like a mean 
flow). 

2.2.2 Axisymmetry induced by an axial vector 

We first consider the case of an axisymmetry characterized by an axial vector n. 
An axisymmetric tensor of rank 4 can only be a function of this vector n and of the 
unit tensor Its most general form, compatible with the symmetry between the 
first two indices, reads: 
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Tliki = bif>ijf>ki + h (f>ik 6 ]i 4 - 6 ,i 6 jk ) 

+ h 8 t} n k rn + b 4 ninjS kl + b 5 ( 6 , k n } ni + 6 jk nini) (2.10) 

+ h {biinjnt t + Sjin,n k ) + b7n t n } n k ni . 

Imposing the constraints ( 2 . 5 ) and defining 62 = — ci, b$ = —ci and 67 = — C3 lead 
to the following expressions: 

61 = (6cj — 4 c 2 n 2 — c 3 n 4 )/ 9 , 

h = b 4 = (4c 2 + c 3 n 2 )/ 3 , (2.11) 

65 = ~ c 2 ■ 

If the constraints ( 2 . 6 ) are imposed on pijkt, only two parameters are different 
from zero and are opposite (65 = -b 6 ). Thus, by introducing 65 = c 4 in p, the 
subgrid-scale stress reads: 


( 2 . 12 ) 

- c 3 - 2c 4 (r.n^ + rurj) , 

where s, = S,*™* and r t = The effect on the resolved energy balance of the 

first three terms is fully determined by the sign of the parameters ci, c 2 , and c 3 . 
Indeed, these terms correspond to dissipation (resp. creation) of resolved energy if 
and only if C\ , c 2 , and c 3 are positive (resp. negative). On the contrary, the sign of 
the term proportional to c 4 in the resolved energy balance depends simultaneously 
on the sign of c 4 and on the flow through the factor 

r*jSij = -ci |S| 2 - 4 c?s 2 - c 3 (s*n*) 2 - 4c 4 sjtrjt . (2.13) 

If the anisotropy is weak (n is relatively small), only terms up to n 2 must be retained; 
since s,,rj = 0(n), the term proportional to c 3 can be neglected in this case. 

2.2.3 Axisymmetry induced by a pseudovector 

We now consider that the anisotropy direction is represented by a pseudovector 
p. The most general axisymmetric tensor of rank 4 will be a function of the vector 
Pi, the unit tensor 6 ij, and the Levi-Civita tensor The situation is thus more 
complicated and more parameters need to be introduced. The notations will be 
simplified by introducing the antisymmetric tensor V l3 = tijkPk so that the most 
general tensor compatible with the symmetry between the first two indices reads: 

Tijki = d\ fiijbki + d 2 (SikSji + S t iSjk) + d$6ijVki 4- d 4 (SikVji + 6jkV u) 

+ d 5 {6,iVj k + SjiVik) + de (eiuPj + tjkiPi) + d 7 6 l jp k pi + d s piPjS k i 

+ d 9 ( 6 tk pjpi + 6 jk p,pi ) + dio ( buPjPk + 6jip,p k ) (2.14a) 

+ d n ( V ik Vji + VuVjk) + du ( V, k p } pi + V ]k p,pi ) 

+ d \ 3 (Vupjpk + Vjipipi t ) + di 4 p,pjV kl + d 15 p,pjp k p, . 


= — 2 ci 5 ij — 2c 2 


yriiSj 


+ SiTij - -6 tJ s k n k 


) 
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We will not discuss the complete tensor T" with 15 independent parameters. Let 
us assume that the anisotropy is weak enough to keep only terms proportional to 
the vector p,. In this case, T n reduces to 

^zjkl ^ d\S t jSkl ~h d 2 

+ dzSjjVki + d 4 (6 ik V }l + 6jkVu) (2.146) 

+ 4 (SuVjk + SjiVik) + de (tzkiPj + tjkiPi) . 

Imposing the constraints (2.5) and defining d 2 = — d$ — —t 2 and c/g = — e 3 lead 
to the following expression: 


d\ — /3 , 

d 3 =2e 3 , (2.15) 

d 4 = — 2e 3 — C '2 , 

while imposing the constraints (2.6) with the new definition c/ 5 = — e 4 , and dg = — e 5 
leads to 

d\ — c /2 0 , 

d 3 = (2e 5 — 4e 4 )/3 , (2.16) 

d 4 — c 4 . 

The subgrid scale stress thus reads: 

r*j = -2e\ S 2J - 2(e 2 + c 3 ) (S^V^- + 

/- _ 2 - \ (2.17a) 

+ 2(e 4 + e 5 ) ( + RjkVik - -SijRkiVki) • 

Although the total eddy viscosity contains 5 parameters, only three of them appear 
independently in the expression for t*j. Let us note that the expression of r t * can 
be simplified by using the resolved vorticity: 

t* = -2ei5 i; - 2(e 2 + e 3 ) (SikVjk + SjkVik) 

/ _ 2 _ \ (2.176) 

+ 2(e 4 + c 5 ) I WiPj + wjp t - -Sijwkpk J . 

It is interesting to note that the anisotropic corrections appear at first order in the 
anisotropy direction p t . Thus, we conclude that a pseudovector anisotropy (like a 
rotation or a magnetic field) should affect the eddy viscosity more rapidly than an 
axial vector anisotropy (like a grid or a flow anisotropy). 

2.3 Onsager reciprocal symmetries 

Strictly speaking, the Onsager reciprocal symmetries do not apply to turbulence. 
Indeed, they have been derived for describing the irreversible return to equilibrium 
in macroscopic system, and they strongly rely on the microscopic reversibility of 
particle motions as well as on the linearity of the transport laws. However, in an 


254 


D. Carati & W . Cabot 


attempt to simplify the eddy viscosity picture as much as possible, it is tempting 
to assume the existence of such relations for the tensors v and ji. We will not try 
to justify further the use of such relations and present the form of the eddy viscos- 
ity tensors fulfilling these relations as a approximate simplification. The Onsager 
reciprocal relation will imply the following additional relations: 

Vijki( n) = t/uiji n), 

Pijkiin) ” Vkiiji n) , ^ ^ 

p) = V k Uj{- p), 

P) — ftkliji P ) • 

When applied to the previous results, these relations imply c 4 = 0 and 65 = — e 4 . 
Thus, they strongly simplify the tensor fi tJ ki but they do not affect the tensor 1 
in the case presented here. 

3* Anisotropic eddy viscosity and dynamic model 

It has been mentioned in the introduction that the use of the dynamic proce- 
dure gives a direct access to a multiple-parameter eddy viscosity. In this section 
we present the dynamic derivation of the eddy viscosity tensor in the simplest 
anisotropic geometry: the weak axisymmetric anisotropy induced by an axial vec- 
tor. Moreover, the problem is further simplified by assuming the existence of On- 
sager symmetries for the tensor v X jki and fi x jkh In that case, we have shown in the 
previous section that the subgrid stress tensor reduces to 

T*j = - 4 c 2 n 2 s[ l J , (3.1) 

where 

sf, = 5^2 i n iSjini + njSuni) - — Sjt/run/^ . (3.2) 

With the new tensorf S^j = Su - 4 and the parameters v\ — c\ + 2n 2 C2 and 
v 2 = Cj, the subgrid stress tensor may be rewritten as 

T : i = -2v i s'! j -2v 2 st J . (3.3) 

With this formulation, the resolved energy dissipation reads e = —T*jSij = U\R\ + 
z^2 R2 1 where 


~ S'jSij - n 2 - ® 1 (3.4a) 

0 

R, = £ S.,f ; = (f,> -§!' ) = 1|5| J - ^ > 0. (3.46) 

IJ ij ' ' 

f This notation should not lead to the conclusion that sf~j and Sjj are orthogonal. It is easy to 
show that Ylij ® m general. 
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The last inequality is a direct consequence of the Cauchy-Schwartz inequality : 

» 2 = BE 5 '*"*)’ s EE (&) (£»?) = »w/ 2 - < 35 > 

l k l k \ l / 

Sufficient conditions for having a positive resolved energy dissipation are thus v\ > 0 
and V 2 > 0- In order to devise the simplest dynamic procedure, we suppose that 
both V\ and u 2 scale following the Kolmogorov law: 


i/i = -CjA 4 / 3 , (3.6a) 

t/ 2 = -C 2 A 4/3 . (3.66) 

The choice for the length scale A in C\ and C 2 (which are not dimensionless) is 
unimportant because the dynamic model will take care of the amplitudes. Only the 
power 4/3 is important. With these definitions, the model becomes 

T tj ~ C\pij 4- CV/i; , (3.7) 

where 

Pij = -2A 4/3 sj l j , (3.8a) 

r h] = -2A 4 ' 3 f£ ■ (3.86) 

Assuming a volume- averaged version of the dynamic model, the error with respect 
to the Germano identity is given by (Germano, 1992; Ghosal et a/., 1995; Lilly, 
1992): 


Eij(C\,C 2 ) = Lij -j- C\Mij + C 2 Nij , (3.9) 

where 

Mij = -2 A 4/3 (1 - a 4/3 )^|' , 

(3-10) 

Nij = -2 A 4 / 3 (1 - a *t 3 )S ijt 

where a is the ratio between test and grid filters. By minimizing Efj , we have the 
two coupled equations: 


{Lij Mu) + C\ {M^ Mij) + C 2 {N tJ M tJ ) = 0 , 
(LiiNij) + C\{MijNij) + C 2 {N l3 N lJ ) - 0. 


Since M l j is not aligned with Nij, these two equations are not linearly proportional 
and they may be used for determining both C\ and C 2 . 
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FIGURE 1. Comparison of the grid anisotropy described in terms of the grid 
spacing in each direction vs. the flow anisotropy described in terms of the square 
root of the diagonal components of the Leonard tensor. Clearly, the flow anisotropy 
is only important close to the wall in the streamwise direction. The grid anisotropy 
is also more important close to the wall, but mainly in the walbnormal direction. 


a) Dx+: ; Dy+: ; Dz+: . b) Square root of Lxx: ; square root 

of Lyy: ; square root of Lzz: . 


4. Application to the channel flow 

The dynamic formulation presented in the previous section has been implemented 
for channel flow with a friction Reynolds number of 1030 (cf. Cabot, 1994). The 
weak axisymmetric anisotropy is probably a very rough approximation for the chan- 
nel geometry, so that the results presented here must be regarded as very prelimi- 
nary tests. Moreover, it is not clear in channel flow which direction is the dominant 
anisotropic one (see Fig. 1). Indeed, channel flow is characterized by two anisotropic 
directions: the streamwise and the wall-normal directions. Both choices for n have 
been tested. 

The rms values of streamwise velocity component (u') are presented in Fig. 2. 
The rms values of the wall-normal and spanwise velocity components (i/ and w') 
seem to be insensitive to the model and are not shown here. It appears that the 
results from the isotropic model and the anisotropic model based on the wall-normal 
direction are almost indistinguishable. The results for the anisotropic model based 
on the streamwise direction seem to be better close to the wall. This could indicate 
that the flow anisotropy has more influence than the grid anisotropy. However, no 
definitive conclusion can be made since the model based on the streamwise direction 
does not perform well in the core region. Also, preliminary results indicate a long- 
time lack of stability for this latter model. 

Finally, we present the results for the two eddy viscosity coefficients (iq and iq) 
in Fig. 3. For both models, the condition of positive dissipation (iq > 0 and iq > 0) 
are mostly well satisfied. It is not yet known if the weakly negative values of iq in 
the model based on the streamwise direction are responsible for its lack of stability. 
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FIGURE 2. Comparison of profiles for the rms of streamwise velocity (id) between 
experimental data, the isotropic dynamic model and two versions of the anisotropic 
dynamic model based on the streamwise and wall-normal anisotropy directions. 

Experimental data (Hussain & Reynolds, 1970):« ; isotropic LES: ; anisotropic 

LES (n=str.dir.): ; anisotropic LES (n=wall norm, dir.): . 



FIGURE 3. Comparison of eddy viscosity coefficients in three different models. 

Isotropic model: v — t ; anisotropic model (streamwise): v\ , v 2 ■; 

anisotropic model (wall normal): u\ , v 2 ■ 

5. Discussion 

The use of an anisotropic eddy viscosity model has been shown to complicate 
dramatically the relation between the subgrid stress tensor and the resolved velocity 
derivatives. In particular, in the fully anisotropic geometry, 40 independent effective 
transport coefficients must be introduced. However, when some approximations are 
used, it is possible to simplify the problem drastically. As an example, we have 
tested the weakly anisotropic axisymmetrical geometry. In that case, the eddy 
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viscosity tensor reduces to a two parameter quantity. A dynamic procedure has 
been proposed for this problem and some tests have been made in channel flow. 

These numerical tests have clearly shown that the determination of the anisotropy 
direction remains an important issue in the simplified anisotropic model presented 
in §3. Indeed, even when the flow is fully anisotropic, the model discussed in §3 may 
be regarded as the first tensorial invariant correction to the isotropic eddy viscosity. 
The use of this model could then be seen as the result of a “local axisymmetric 
assumption” which should be at least as robust as the local isotropic assumption. 
However, in that case it is probably crucial to chose the vector n in an appropriate 
way. It is also possible that the vector n varies with space. An interesting exten- 
sion to this work would be the derivation of a dynamic procedure giving explicit 
expressions not only for the eddy viscosity amplitudes but also for the vector n. 

At this point the simplest test for anisotropic models would be the homogeneous 
rotating turbulence. In that case, the anisotropy direction is clearly determined and 
is given in terms of the rotation pseudovector. 
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